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Dynamics  and  Control  of  a  Minimally  Actuated  Biomimetic 

Vehicle: 

Part  I  -  Aerodynamic  Model 

David  B.  Doman  * 

Michael  W.  Oppenheimer  t 
David  O.  Sigthorsson  + 


An  aerodynamic  model  for  the  forces  and  moments  acting  on  a  minimally 
actuated  flapping  wing  micro  air  vehicle  (FWMAV)  are  derived  from  blade 
element  theory.  The  proposed  vehicle  is  similar  to  the  Harvard  RoboFly 
that  accomplished  the  first  takeoff  of  an  insect  scale  flapping  wing  air¬ 
craft,  except  that  it  is  equipped  with  independently  actuated  wings  and 
the  vehicle  center-of-gravity  can  be  manipulated  for  control  purposes.  Us¬ 
ing  a  blade  element-based  approach,  both  instantaneous  and  cycle-averaged 
forces  and  moments  are  computed  for  a  specific  type  of  wing  beat  motion 
that  enables  nearly  decoupled,  multi-degree-of-freedom  control  of  the  air¬ 
craft.  The  wing  positions  are  controlled  using  oscillators  whose  frequencies 
change  once  per  wing  beat  cycle.  A  new  technique  is  introduced,  called 
Split-Cycle  Constant-Period  Frequency  Modulation,  that  has  the  desirable 
property  of  providing  a  high  level  of  control  input  decoupling  for  vehicles 
without  active  angle-of-attack  control.  Like  the  RoboFly,  the  wing  angle-of- 
attack  variation  is  passive  by  design,  and  is  a  function  of  the  instantaneous 
angular  velocity  of  the  wing  in  the  stroke  plane.  A  control-oriented  dy¬ 
namic  model  of  the  vehicle  is  derived,  which  is  based  on  a  cycle-averaged 
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representation  of  the  forces  and  moments.  Control  derivatives  are  calcu¬ 
lated  and  a  cycle-averaged  control  law  is  designed  that  provides  control  over 
5  degrees-of-freedom  of  the  vehicle. 

I.  Introduction 

The  first  takeoff  of  an  insect-scale  flapping  wing  micro  air  vehicle  was  achieved  by  an 
aircraft  called  RoboFly  that  was  developed  at  Harvard  University  by  Wood  et.al.1  A  key 
feature  that  led  to  the  successful  first  flight  is  that  the  vehicle  was  minimally  actuated. 
RoboFly  uses  a  single  bimorph  piezoelectric  actuator  to  impart  symmetric  stroke  plane 
motion  to  two  wings.  The  wing  angle-of-inclination  with  respect  to  the  stroke  plane  is 
passively  regulated  by  wing  flexure  joints  resulting  in  planform  motion  that  is  similar  to 
that  of  a  dipterian  insect.  This  passive  rotation  of  the  planform  about  the  spar  eliminates 
the  need  for  added  weight  and  complexity  associated  with  actuators  that  actively  rotate  the 
planform  about  the  spar.  The  first  flight  at  Harvard  resulted  in  unregulated  flight  up  a  pair 
of  wires  that  constrained  the  vehicle  motion  to  vertical  translation. 

In  the  present  paper,  a  vehicle  concept  and  a  control  strategy,  that  would  enable  con¬ 
trolled  six-degree-of-freedom  flight  of  the  fuselage  of  an  aircraft  similar  to  the  RoboFly 
without  the  need  for  artificial  motion  constraints,  are  developed.  A  diagram  of  the  proposed 
vehicle  is  shown  in  Figure  1. 
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Figure  1.  General  assembly  of  a  minimally  actuated  flapping  wing  micro  air  vehicle. 


The  differences  between  the  proposed  aircraft  and  the  Harvard  RoboFly  are  that  it  would 
be  equipped  with  independently  actuated  wings  and  an  actuated  bob-weight  that  would  allow 
the  vehicle  center-of-gravity  to  be  manipulated  for  control  purposes.  The  wings  rotate  about 
hinges  at  the  wing  root.  A  linkage  based  transmission  translates  the  tangential  motion  of 
the  tip  of  a  bimorph  piezoelectric  actuator  into  rotational  motion  of  the  wings  in  the  stroke 
plane.  The  linkage  elements  are  designed  to  achieve  impedance  matching  between  the  wing 
and  actuator  forces  as  well  as  amplify  the  relatively  small  motion  of  the  tip  of  the  bimorph 
strip  into  large  angular  displacements  of  the  wing  root.  The  planforms  are  connected  to  the 
movable  wing  roots  by  a  limited  hinge  joint  that  provides  for  passive  rotation  of  the  wing. 
This  hinge  allows  the  wing  to  passively  flip-over  as  the  wing  reverses  direction  at  the  end 
of  each  stroke  and  enables  the  chord  to  rotate  about  the  axis  of  the  spar  in  a  manner  that 
approximates  the  wing- twisting  motion  that  has  been  observed  in  dipterian  insects.2  As  the 
wing  rotates  through  the  stroke  plane,  dynamic  pressure  acting  on  the  wing  tends  to  cause  it 
to  feather  into  the  wind;  however,  as  shown  in  Figure  2,  the  spar  and  root  structure  interfere 
at  an  angle  set  by  the  designer  to  prevent  the  wing  from  over-rotating.  This  interference 
causes  the  wing  to  hold  a  constant  angle-of-attack  relative  to  the  stroke  plane  once  a  critical 
dynamic  pressure  is  reached. 
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Downstroke 


End  of  Stroke 


Upstroke 


Figure  2.  Detail  of  passive  wing  rotation  joint. 


The  bimorph  piezoelectic  actuators  and  the  carbon  fiber  substrate  to  which  they  are 
mounted  are  cantilevered  to  the  fuselage.  In  the  case  of  the  wings,  the  actuator  assembly 
and  linkage,  together  with  the  wing,  form  a  spring-mass-damper  system  that  has  a  known 
resonant  frequency.  The  actuator  assembly  for  the  bob-weight  is  simply  used  to  control 
the  mean  pitch  attitude  of  the  fuselage.  In  the  single  actuator  Harvard  experiments,  the 
dynamic  system  was  driven  at  resonance3  for  maximum  energy  efficiency  to  achieve  flight. 
In  the  present  problem,  each  wing-linkage-actuator  system  is  nominally  driven  at  a  hover 
frequency  that  is  defined  as  the  frequency  at  which  the  cycle-averaged  lift  is  equal  to  the 
weight  of  the  aircraft.  The  term  hover  is  to  be  interpreted  in  a  time-averaged  sense  because 
as  we  will  show,  the  aircraft  is  in  constant  motion  due  to  the  periodic  nature  of  the  forces 
and  moments  produced  by  the  flapping  wings.  The  closest  approximation  of  hover  that  can 
be  achieved  is  that  of  a  high  frequency,  low  amplitude  limit  cycle  about  a  mean  position.4 

Using  a  blade  element  based  aerodynamic  model  and  time  averaging,  it  will  be  shown 
that  5  independent  degrees  of  freedom  can  be  controlled  and  that  the  6th  degree-of- freedom 
can  be  manipulated  by  coordinated  motion  of  a  subset  of  the  5  degrees  of  freedom.  The  wing 
angle-of-attack  variation  is  passive  and  is  a  function  of  the  instantaneous  angular  velocity  of 
the  wing  in  the  stroke  plane.  The  control  laws  are  therefore  designed  to  make  use  of  three 
actuators,  two  of  which  control  the  angular  position  of  the  wing  in  the  stroke  plane  and  one 
that  moves  a  bob-weight  that  manipulates  the  vehicle  center-of-gravity.  A  technique  called 
Split-Cycle  Const  ant- Period  Frequency  Modulation  is  proposed  that  allows  each  wing  to 
generate  non-zero  cycle-averaged  rolling  and  yawing  moments  as  well  as  translational  forces 
that  can  accelerate  the  vehicle  forwards  and  backwards  while  maintaining  a  wingbeat  forcing 
function  that  keeps  the  wing-actuator  dynamics  operating  near  the  hover  frequency.  The 
split-cycle  technique  achieves  this  objective  by  varying  the  frequency  of  the  oscillators  that 
drive  the  wing  position  actuators,  such  that  the  dynamic  pressure  acting  on  the  wing  during 
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the  forestroke  is  different  from  that  acting  on  the  aftstroke.  Pitching  moment  and  high  speed 
translation  are  controlled  by  varying  the  vehicle  center- of- gravity. 

A.  Preliminary  Assumptions 

The  purpose  of  this  investigation  is  to  derive  an  idealized  model  that  will  allow  one  to  study 
the  dynamic  behavior  of  the  aircraft  under  consideration  and  to  formulate  and  test  control 
strategies  using  non-real-time  simulations.  This  objective  limits  the  fidelity  of  the  aerody¬ 
namic  and  structural  models  that  can  be  used.  The  use  of  computational  fluid  dynamics  and 
finite  element  structural  codes  would  be  computationally  intractable  for  a  simulation  of  this 
type  and  would  likely  impede  physical  insights  that  would  allow  one  to  formulate  a  control 
law  and  observe  key  vehicle  design  features.  In  light  of  the  above  objective,  the  following 
assumptions  are  made: 

•  There  are  no  aerodynamic  interactions  between  the  left  and  right  wings. 

•  There  are  no  aerodynamic  interactions  between  the  wings  and  the  fuselage. 

•  The  2D  sectional  aerodynamic  coefficients  are  known  and  constant  throughout  each 
stroke.  Three  dimensional  and  unsteady  effects  from  leading  edge  vortex  and  wake 
phenomena  are  not  present. 

•  The  passive  wing  rotation  joint  is  on  a  limit  when  the  wing  angular  velocity  in  the 
stroke  plane  is  non-zero. 

•  The  air  mass  surrounding  the  vehicle  is  quiescent. 

•  Aerodynamic  forces  and  moments  are  the  sole  result  of  wing  motion,  i.e.,  only  flight 
conditions  in  the  neighborhood  of  hover  are  considered. 

•  The  bandwidth  of  the  piezoelectric  actuators  exceeds  the  hover  frequency. 

II.  Coordinate  Frame  Definitions 

The  dynamic  analysis  of  the  vehicle  requires  the  use  of  several  coordinate  frames.  This 
is  because  the  signs  of  aerodynamic  forces  and  moments  differ  from  the  right  wing  to  the 
left.  Furthermore,  the  signs  of  the  forces  and  moments  are  dependent  upon  the  direction 
of  the  wing  stroke,  i.e.,  upstroke  or  downstroke.  The  aerodynamic  forces  and  moments 
can  be  conveniently  written  in  certain  intermediate  frames;  however,  it  will  be  necessary  to 
transform  these  parameters  into  a  body-fixed  and  ultimately  an  inertial  coordinate  frame 
in  order  to  write  the  equations  of  motion  of  the  fuselage.  Definitions  of  the  numerous 
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coordinate  frames  and  rotation  matrices  required  in  this  analysis  are  provided  in  Tables  1 
and  2,  respectively.  The  coordinate  frame  and  angle  definitions  are  illustrated  in  Figure  3. 


Figure  3.  Relationship  between,  body,  root,  spar,  left  and  right  wing  planform  axis  systems 
on  upstroke. 


A  system  of  notation  has  been  selected  where  superscripts  “cancel”  subscripts  when 
multiplying  rotation  matrices.  For  vector  quantities,  superscripts  identify  the  coordinate 
frame  in  which  the  quantities  are  expressed,  while  subscripts  identify  the  quantity. 
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Frame 

Symbol 

Feature 

Location/Sense 

Body 

(Fixed  to  Fuselage) 

B 

Origin 
xB  axis 
yB  axis 
zB  axis 

Nominal  Center  of  Gravity 
(+)  Anterior  Normal  to  Stroke  Plane 
(+)  Right  Hand  Side 
(+)  Ventral  Side 

Inertial 

(Fixed  to  Flat  Earth) 

I 

Origin 
xi  axis 
yi  axis 
zj  axis 

Initially  Coincident  with  Body  Center-of-Gravity 
(+)  Initially  Aligned  with  xr  axis 
(+)  Initially  Aligned  with  yB  axis 
(+)  Initially  Aligned  with  zr  axis 

RW  Root 

(Fixed  to  Fuselage) 

RWR 

Origin 
xrwr  axis 
Vrwr  axis 
zrwr  axis 

RW  Root  Hinge  Point 
(+)  in  Xr  dir. 

(+)  in  yB  dir. 

(+)  in  zB  dir. 

LW  Root 

(Fixed  to  Fuselage) 

LWR 

Origin 
xlwr  axis 
Ulwr  axis 
zlwr  axis 

LW  Root  Hinge  Point 
(+)  in  Xr  dir. 

(+)  in  yB  dir. 

(+)  in  zB  dir. 

RW  Spar 

(Rotates  with  RW  Spar) 

RWS 

Origin 

XRWS  axis 
Urws  axis 
zrws  axis 

RW  Root  Hinge  Point 

(+)  Towards  ventral  side  at  (j)Rw  =  0 

(+)  Proximal  to  distal,  coincident  c  RW  spar 

(+)  Towards  posterior  of  fuselage 

LW  Spar 

(Rotates  with  LW  Spar) 

LWS 

Origin 
xlws  axis 
Vlws  axis 
zlws  axis 

LW  Root  Hinge  Point 
(+)  Towards  dorsal  side  at  4>lw  =  0 
(+)  Proximal  to  distal,  coincident  c  spar 
(+)  Towards  posterior  of  fuselage 

RW  Planform  Upstroke 
(Aligned  c  RW 

Plane  on  Upstroke) 

RWPU 

Origin 
xrwpu  axis 
Vrwpu  axis 
zrwpu  axis 

RW  Root  Hinge  Point 

(+)  Direction  of  wing  rotation 

(+)  Proximal  to  distal,  coincident  c  RW  spar 

(+)  anterior,  completes  RH  coord,  sys 

RW  Planform  Downstroke 
(Aligned  c  RW 

Plane  on  Downstroke) 

RWPD 

Origin 
xrwpd  axis 
Vrwpd  axis 
zrwpd  axis 

RW  Root  Hinge  Point 

(+)  Direction  of  wing  rotation 

(+)  Proximal  to  distal,  coincident  c  RW  spar 

(+)  posterior,  completes  RH  coord,  sys 

LW  Planform  Upstroke 
(Aligned  c  LW 

Plane  on  Upstroke) 

LWPU 

Origin 
xlwpu  axis 
Ulwpu  axis 
zlwpu  axis 

LW  Root  Hinge  Point 

(+)  Direction  of  wing  rotation 

(+)  Proximal  to  distal,  coincident  c  LW  spar 

(+)  posterior,  completes  RH  coord,  sys 

LW  Planform  Downstroke 
(Aligned  c  LW 

Plane  on  Downstroke) 

LWPD 

Origin 
xlwpd  axis 
Ulwpd  axis 
zlwpd  axis 

LW  Root  Hinge  Point 

(+)  Direction  of  wing  rotation 

(+)  Proximal  to  distal,  coincident  c  LW  spar 

(+)  anterior,  completes  RH  coord,  sys 

Table  1.  Coordinate  Frame  Definitions. 
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From 

To 

Stroke 

Symbol 

Rotation  Matrix 

RW  &  LW  Root 

Body 

Either 

K 

jdB  _  tdB  _  T 

nRWR  ~  rt LWR  ~  1 

RW  Spar 

RW  Root 

Either 

tdRWR 

11 RWS 

0  0  -1 

-  sin  (pnw  cos  <j)RW  0 
_  cos  (pnw  Sin  (pnw  0  _ 

LW  Spar 

LW  Root 

Either 

r>LWR 

11 LWS 

[ 

0  0  -1 

—  sin  (f>LW  —  cos  4>lw  0 

—  cos  (j)LW  sin  <i)LW  0 

1 

RW  Planform  Upstroke 

RW  Spar 

Up 

tdRWS 

11 RWPU 

—  cos  a  0  sin  a 

0  1  0 

—  sin  a  0  —  cos  a 

RW  Planform  Downstroke 

RW  Spar 

Down 

tdRWS 

RWPD 

cos  a  0  sin  a 

0  1  0 

—  sin  a  0  cos  a 

LW  Planform  LIpstroke 

LW  Spar 

Up 

jdLWS 

LWPU 

cos  a  0  sin  a 

0  1  0 

—  sin  a  0  cos  a 

LW  Planform  Downstroke 

LW  Spar 

Down 

jdLWS 

11 LWPD 

—  cos  a  0  sin  a 

0  1  0 

—  sin  a  0  —  cos  a 

Table  2.  Transformation  Matrices. 


III.  Aerodynamic  Model  of  Flapping  Wing  Vehicle  with 

Independently  Actuated  Wings 

A.  Instantaneous  Aerodynamic  Forces  and  Centers  of  Pressure  in  Wing  Plan- 
form  Frames 

The  aerodynamic  forces  are  derived,  using  blade  element  theory,  for  triangular  shaped  wings 
that  have  two  degrees  of  freedom,  namely,  angular  displacement,  </>(£),  about  the  wing  root  in 
the  stroke  plane,  and  angular  displacement  of  the  planform  about  the  passive  rotation  hinge 
joint,  which  is  equivalent  to  wing  angle-of-attack,  a,  in  still  air.  The  triangular  planform 
wing,  shown  in  Figure  4,  is  taken  to  be  a  rigid  flat  plate  whose  elemental  lift  and  drag  at  a 
spanwise  location  ywp  is  given  by: 

dL  =  ^CL(a)<j)2y^VPc(yWp)dyWP  (1) 
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(2) 


dD  =  ^ CD(oi)4>2ywpc(ywp)dywp 

where  c(ywp )  is  the  chord  at  the  spanwise  location  yWp ,  which  is  a  location  on  the  wing 
spar.  Now  compute  the  total  lift  and  drag  on  the  wing  expressed  in  a  wing  planform  fixed 
coordinate  system,  e.g.,  RWPU,  RWPD,  LWPU,  and  LWPD.  Such  a  coordinate  frame  has  an 
origin  at  the  wing  root  hinge  point  and  the  x-y  plane  is  coincident  with  the  wing  planform. 
The  lift  is  computed  by  integrating  the  elemental  lift  over  the  span  according  to 


n  v 


WP 


yWP 


Figure  4.  Blade  element  computation  of  aerodynamic  forces,  moments  and  centers  of  pressure. 


L=  [R  dL  =  ^CL(a)^(t)2IA  (3) 

Jo  ^ 

Similarly,  drag  is  computed  according  to 

D=  [R  dD=%CD(a)j>(t)2IA  (4) 

Jo  ^ 

where  I  a  is  the  area  moment  of  inertia  of  the  planform  about  the  root  and  R  is  the  length  of 
the  wing.  For  convenience  all  of  the  time  invariant  parameters  are  lumped  together  according 
to 

kL  =  P- CL(a)IA 
kD  =  ^ CD(a)IA 
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(5) 


thus,  lift  and  drag  can  be  expressed  as  the  product  of  time  invariant  parameters  and  time 
varying  functions 


L  =  kL(j)(t)2 

(o) 

D  =  kD<j*(t )2 

The  following  expressions  for  lift  and  drag  coefficients,5  which  were  obtained  from  low 
Reynolds  Number  experiments,  are  used  in  this  work 


CD(a)  =  1.92  -  1.55cos(2.04a  -  9.82) 
CL(a)  =  0.225  +  1.58 sin (2. 13a  -  7.2) 


where  a  is  in  degrees.  Note  that  the  only  variable  that  can  be  actively  manipulated  to  control 
the  instantaneous  aerodynamic  forces  is  the  angular  velocity  of  the  wing  0(f)  and  that  the 
forces  are  quadratic  functions  of  this  motion  variable.  The  center  of  pressure  in  the  plane 
of  each  wing  must  be  obtained  in  order  to  compute  moments  on  the  fuselage.  The  center 
of  pressure  location  for  each  wing  in  their  respective  local  wing  planform  coordinate  frame, 
namely  RWPU,  RWPD,  LWPU,  or  LWPD,  are: 


x 


WP 

cp 


y 


WP 

cp 


_ WP 

^ cp 


Jo  VwP2c(ywp)~dywp  — c  / 3b2  +  2bR  +  R2 \ 

f0R  ywpc(ywp)dyWp  5  v  b2  +  bR  +  R2  J 

Jo  y\ypc(ywp)dywp  _  3  /  Id  +  b2 R.  +  bR 1  +  Ri\ 
f0*  ywpc(ywp)dyWp  5  V  3 b2  +  2 bR  +  R2  J 
0 


(8) 

(9) 

(10) 


where  WP  denotes  any  of  the  four  local  wing  planform  coordinate  systems,  b  is  the  distance 
from  the  wing  root  to  the  wing  break  along  the  ywp  axis,  and  c  is  the  maximum  chord  of 
this  wing. 


B.  Expression  of  Aerodynamic  Forces  and  Centers  of  Pressure  in  Body  Frame 

With  the  relationships  between  the  body,  roots,  spars,  upstroke  planform  and  downstroke 
planform  axis  systems  established  in  Tables  1  and  2,  the  instantaneous  values  of  lift  and 
drag  on  each  wing  are  transformed  into  the  body-axis  coordinate  frame.  Recall  that  lift  is 
defined  as  the  component  of  aerodynamic  force  perpendicular  to  the  relative  wind,  while 
drag  is  defined  as  the  component  of  aerodynamic  force  parallel  to  the  relative  wind.  If 
the  air  mass  is  quiescent,  then  the  relative  wind  is  parallel  to  the  stroke  plane,  which  is 
coincident  with  the  x-y  planes  of  the  local  spar-frames.  Therefore,  the  lift  and  drag  forces 
are  conveniently  expressed  in  the  spar  coordinate  frames,  which  can  be  transformed  to  the 
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body  frame  using  the  rotation  matrices  defined  in  Table  2.  For  example,  to  express  the 
aerodynamic  force  vector  associated  with  the  right  wing  on  the  upstroke  in  the  body-frame, 
the  following  transformation  is  applied 


F 


B  _ 
RWU  ~ 


t-)  b  R  rwrt?rws 

rvRWRrvRtVSr  RWU 


(ii) 


The  aerodynamic  forces  associated  with  each  wing  and  stroke  expressed  in  both  the  spar 
and  body  frames  are  summarized  in  Table  3. 


Force 

Local  Spar  Frame 

Body  Frame 

RW  Upstroke 

t?RWS  _ 
r  RWU  — 

Drwu 

0 

—  LrwU_ 

-pi? 

r  RWU  — 

Lrwu 

— Drwu  sin  4>R.wu{t ) 

_  DrWU  cos  <pRwu(t)  _ 

RW  Downstroke 

rpittVS  _ 
r  RWD  — 

—  DrwD 

0 

_  —  LrwD_ 

t?B 

r  RWD  — 

Lrwd 

Drwd  sin  4>rwd  if) 

—Drwd  COS  (f)RW o(t) 

LW  Upstroke 

T71 LWS 
r  LWU  ~ 

— dlwu 

0 

—Llwu_ 

T7 \B 

r  LWU  ~ 

Llwu 

Dlwu  sin  <puwu(t) 
_Dlwu  cos  4>Lwu(t)_ 

LW  Downstroke 

t?LWS  _ 
r  LWD  ~ 

Dlwd 

0 

—Llwd_ 

T pi? 

r  LWD  ~ 

Llwd 

—  Dlwd  sin  ( t>Lwo(t ) 
—Dlwd  cos  4>LWD(t)_ 

Table  3.  Aerodynamic  forces  expressed  in  local  spar  and  body  frames. 


The  center  of  pressure  on  each  wing  is  conveniently  expressed  in  the  local  wing  planform 
frame  associated  with  each  stroke.  To  calculate  the  center  of  pressure  in  the  body  frame 
three  coordinate  frame  rotations  and  one  translation  must  be  performed.  One  rotation  is 
trivial  because  the  root  and  body  systems  are  parallel.  The  translation  is  associated  with 
the  distance  between  the  origin  of  the  body  frame  and  the  origin  of  each  wing  root  frame. 
For  example,  for  the  right  wing  upstroke,  the  center  of  pressure  location  in  the  body  frame 
is  given  by 


B 


„  _  TD  B 

lcPRWU  ~  *xRWRtxRWS*xRWPULcPRWU 


,RWR-r,RWS 


RWPU 


+  Ar 


(12) 


where  Ar^  is  the  position  vector  from  the  origin  of  the  body  axes  coordinate  system  to  the 
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origin  of  the  right  wing  root  coordinate  system,  i.e., 


Arg 

Arf 


A 

=  r 


B 

°RWR 


A 

=  r 


B 

°LWR 


Ax%  f  A  z% 
A xf  —  ^  A  zp 


(13) 


The  width  of  the  vehicle  is  designed  as  w  and  the  origin  of  the  body  axis  coordinate  system 
is  assumed  to  be  located  at  the  midpoint  of  the  fuselage  in  the  y-body  axis  direction.  The 
centers  of  pressure  associated  with  each  wing  and  stroke,  expressed  in  the  body  frame,  are 
summarized  in  Table  4. 


CP  Location 

Body  Frame  Expression 

RW  Upstroke 

B 

1  CPRWU 

x^pp  sin  a  +  Axp 

x%p  sin  (j)RW  cos  a  +  y^p  cos  (j)RW  +  ^ 
-X^pp  cos  (j)RW  COS  a  +  y™p  sin  <pRW  +  A zp 

RW  Downstroke 

Tcprwd  = 

x^p  sin  a  +  Axp 

-x^p  sin  (pRW  cos  a  +  y™p  cos  <j>RW  +  f 
_x^p  cos  (paw  cos  a  +  yWp  sin  <j>RW  +  A z%_ 

LW  Upstroke 

Tcplwu  = 

x^pp  sin  a  +  Axf 

—x^pp  sin  (pLW  cos  a  —  y^p  cos  4>lw  —  f 
rx%P  cos  0LW  cos  a  +  VcpP  sin  & LW  +  ^zl  . 

LW  Downstroke 

rcp  LWD  = 

x^pp  sin  a  +  Axp 

x%p  sin  (f>Lw  cos  a  —  y^p  cos  4>lw~  f 
x^p  cos  4>lw  cos  a  +  y^pp  sin  (j)LW  +  Azf 

Table  4.  Centers  of  pressure  expressed  in  body  frame. 


IV.  Aerodynamic  Moments  in  Body  Frame 


The  expressions  for  the  aerodynamic  moments  associated  with  each  wing  and  stroke  are 
given  by 


^RWU  =  rcPRWU  X  ^RWU 
M-RWD  =  rcp BWD  x  ^rwd 
M-LWU  —  rcpLWU  X  ^LWU 
^-LWD  =  rcpLWD  X  ^LWD 


(14) 
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Carrying  out  the  cross  product  operations  and  substituting  the  values  from  Tables  3  and  4 
into  the  expressions  in  Equation  14  yields 


M 


B 

RWU 


Drwu  \y% F  +  f  cos  <pRW  +  A zp  sin  (j)Rw\ 

{Lrwu  \y^pP  sin  4>rw  +  ^zr]  ~  DrwuAxp  cos  <pRW  —  . . . 

[Lrwu  cos  a  +  Drwu  sin  a ]  x cos  (j)RW  } 
{—DrwuAxx  sin  (j)RW  —  Lrwu  [j  +  y^pp  cos  4>rw\  —  ■  ■  ■ 
[Lrwu  cos  a  +  DRWU  sin  a\  x^p  sin  (pmv  } 


(15) 


M 


B  _ 
RWD  ~ 


M 


B  _ 
LWU  — 


M 


B  _ 
LWD  ~ 


—Drwd  \yWp  +  f  cos  (j)Rw  +  A Zr  sin  <pRw\ 

{Lrwd  \yWp  sin  (p rw  +  A Zr\  +  DrwdAxr  cos  <pRw  +  •  •  • 
[Lrwd  cos  a  +  Drwd  sin  a ]  x^pp  cos  (j)RW } 
{DRwdAxp  sin  (f>mv  —  LRwd  [f  +  y^P  cos  (Pmv]  +  •  •  • 
[Lrwd  cos  a  +  DRWd  sin  a]  x]pp  sin  (pRW} 

Dlwu  [~y%P  -  f  cos  (pLw  -  A zp  sin  (pLW\ 

{ Llwu  [y^P  sin  (pLw  +  A^f  ]  —  DLWu^l  cos  4>lw  —  ■  ■  ■ 
[Llwu  cos  a  +  DLWu  sin  a]  x™pp  cos  (pLw } 

{ DLWu Axf  sin  <pLw  +  LLwu  [f  +  y^pP  cos  4>lw]  +  •  •  • 
[Llwu  cos  a  +  DLwu  sin  a]  x^pp  sin  4>lw} 

Dlwd  [y%P  +  f  cos  (pLw  +  A zp  sin  cpLw\ 

{Llwd  [y^p  sin  <P>lw  +  A^]  +  DlwdAxp  cos  <P>lw  +  •  •  • 
[Llwd  cos  a  +  Dlwd  sin  a]  x^p  cos  (pLw } 

{— DlwdAxp  sin  cpLw  +  LLwd  [f  +  y^P  cos  4>lw\  —  ■  ■  ■ 
[Llwd  cos  a  +  DLwd  sin  a]  x^p  sin  cpLw } 


(16) 


(17) 


(18) 


Equations  15-18  provide  the  expressions  for  the  instantaneous  aerodynamic  moments  gener¬ 
ated  by  each  wing  at  any  point  in  a  wing  beat  cycle  in  the  body-axis  coordinate  frame. 


V.  Split-Cycle  Constant-Period  Frequency  Modulation  Control 

Strategy 

Because  the  vehicle  is  designed  to  mimic  dipterian  insect  flight,  cp  must  be  a  time  varying 
function  that  is  equal  to  zero  at  the  extreme  limits  of  wing  position.  It  is  assumed  that  the 
wing  position  in  the  stroke  plane  can  be  controlled  directly  via  a  high  bandwidth  piezoelectric 
bimorph  actuator.  This  assumption  approximates  the  physics  of  applying  a  voltage  to  a 
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bimorph  piezoelectric  actuator  that  imparts  rotational  motion  to  a  wing  via  a  linkage  that 
consists  of  rigid  elements  connected  by  pin  joints.  The  forcing  function  that  drives  the 
rotation  of  each  wing  is  selected  as 

(f)  =  cos  cut  (19) 

As  defined  in  Figure  3,  the  maximum  value  of  this  angle  occurs  when  either  wing  reaches  its 
maximum  displacement  towards  the  ventral  side  of  the  vehicle.  Assuming  that  the  frequency 
of  the  oscillating  wing  is  held  constant  over  the  segment  of  interest  for  each  wingbeat  cycle, 
the  angular  velocity  of  the  wing  is  given  by 


(p  =  —  u  sin  ut  (20) 

Note  that  the  units  of  </>  are  in  radians  and  that  the  amplitude  of  the  wing  rotation  in  the 
stroke  plane  is  taken  to  be  ±1  rad,  which  closely  approximates  the  rotation  limits  used  in 
the  Harvard  RoboFly.  The  frequency  of  the  oscillator  that  drives  the  actuator  is  selected 
as  the  control  input  variable  to  enable  multi-degree-of- freedom  flight,  and  this  frequency  is 
allowed  to  change  at  key  points  in  the  wingbeat  cycle,  namely  at  the  ends  of  the  stroke 
when  (j)  =  l,(p  =  0  and  (j)  —  —  l,<j)  —  0.  In  a  companion  paper  concerning  1  degree-of- 
freedom  (DOF)  flight,4  it  was  shown  that,  if  the  wingbeat  frequency  is  held  constant  over 
each  cycle,  then  the  cycle-averaged  longitudinal  force  along  the  x-body  axis  was  always  finite 
and  positive,  while  the  cycle-averaged  z-body  force  was  always  zero.  Such  conditions  made 
it  possible  to  achieve  a  practical  condition  that  approximates  hover4  and  allows  altitude 
tracking;  however,  in  order  to  enable  a  vehicle  to  achieve  multi-DOF  flight,  an  additional 
control  input  variable  must  be  introduced.  In  the  interest  of  maintaining  the  simplicity  of  the 
actuation  system,  the  controllability  of  a  vehicle  equipped  only  with  independently  actuated 
wings  with  passive  rotation  about  the  chord  is  examined.  In  order  to  achieve  a  non-zero 
time  averaged  body  force  in  the  z-direction,  subject  to  the  imposed  actuation  constraint, 
the  blade  element  velocities  must  be  asymmetric  over  the  wing  beat  cycle.  Although  there 
are  many  ways  to  accomplish  this,  the  waveforms  that  define  the  motion  of  the  wing  in  the 
stroke  plane  are  parameterized  such  that  the  frequencies  of  the  upstroke  and  downstroke 
can  be  modified  while  maintaining  the  original  period.  This  technique  is  called  Split-Cycle 
Constant-Period  Frequency  Modulation.  Each  wing  will  be  driven  by  a  position  command 
of  the  following  form 


=  cos  [(cu  -  S)t], 
=  cos  [(cu  +  a)t  +  £], 


0  <  t  <  * 

—  (ui-5) 

(21) 

7  ““77  <t<  — 

(uj—o)  —  U) 

(22) 
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where  the  (f>u(t)  and  4>D(t)  represent  the  angular  position  command  associated  with  a  given 
wing  on  the  upstroke  and  downstroke,  respectively,  6  is  the  split-cycle  parameter, 


_a  hen 
c o  —  2  5 


and  the  phase  shift  on  downstroke  is  given  by 


—2n5 
uj  —  25 


(23) 


(24) 


The  phase  shift  enforces  the  compatibility  condition  where  the  frequency  of  the  waveform 
changes,  i.e., 

(cu  +  a)t\t=-2L_  +  f  =  7T  (25) 

Note  that  the  split-cycle  parameter,  5,  and  the  fundamental  frequency,  tu,  are  control  input 
variables  and  that  cr  is  defined  such  that  the  wing  beat  cycle  is  completed  at  exactly  the 
same  time  as  it  would  be  if  the  wing  beat  were  temporally  symmetric,  i.e.,  5  =  a  =  0.  The 
angular  velocity  of  each  wing  over  each  split-cycle  interval  is  given  by 


4>u(t)  =  —  (<n  —  5)  sin  [(cn  —  5)t], 
<j>D(t )  =  —  (<n  +  a)  sin  [(tn  +  a)t  +  £], 


0  <  t  <  ,  * 

—  (u—6) 

(26) 

■7-^  <t<  — 

(LJ  —  d)  —  LO 

(27) 

VI.  Cycle- Averaging 

The  control  strategies  proposed  in  this  paper  are  based  on  the  assumption  that  the  band¬ 
width  of  the  fuselage  controller  is  much  less  than  the  trim  flapping  frequency  required  for 
hover.  If  a  non-oscillatory  control  force  and  inertial  measurements  were  available  for  feed¬ 
back,  one  could  use  conventional  control  design  techniques  to  synthesize  feedback  control  laws 
that  produce  favorable  closed  loop  responses.  Here,  time  varying  high  frequency  oscillatory 
control  inputs  must  be  used;  therefore,  the  relationship  between  the  cycle-averaged  forces 
and  moments  and  the  control  input  parameters,  namely  u>lw,  5lw,  w rw  and  8rwi  which  are 
the  fundamental  frequency  and  split-cycle  parameters  for  the  left  and  right  wings,  respec¬ 
tively,  are  computed.  Feedback  control  laws  based  on  cycle-averaged  forces  and  moments 
will  allow  a  vehicle  to  track  desired  angular  and  spatial  positions  in  a  mean  sense;  however, 
because  of  the  true  periodic  nature  of  the  aerodynamic  forces,  the  vehicle  will  exhibit  limit 
cycle  behavior  in  a  neighborhood  about  the  mean  position.  The  objective  of  this  work  is  to 
develop  controllers  that  drive  the  vehicle  to  a  desired  mean  position  and  orientation. 

The  general  method  used  to  compute  the  cycle-averaged  aerodynamic  forces  and  their 
sensitivities  to  variations  in  the  control  input  parameters  will  now  be  described.  Define 
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G(t)  to  be  a  generalized  force  aligned  with  a  principal  body-axis  direction  that  may  repre¬ 
sent  either  a  force  or  a  moment.  In  order  to  compute  the  cycle-averaged  generalized  force 
associated  with  each  wing,  an  integral  of  the  following  form  is  evaluated 


G 


to 

2n 


(28) 


Since  the  introduction  of  the  split-cycle  parameter  changes  the  frequency  of  the  cosine  wave 
when  (to  —  S)t  =  i r  in  each  cycle,  it  is  convenient  to  split  the  integral  as  follows 


G 


to 

2n 


G(<f>u(t))dt  + 


G(4>D(t))dt 


(29) 


The  control  derivatives  associated  with  a  generalized  time-averaged  force  are  calculated  by 
computing  the  partial  derivative  with  respect  to  each  control  input  variable  about  the  hover 
frequency  toQ. 


dG 

Gauji  ~ 

dtoi 

dG 

GSl  = 

dSi 

dG 

Gau)T  — 

dtOr 

dG 

Gsr  = 

d5r 

uJi=uo,8i=0 


lji=ljo,8i=0 


LUf  — LOq  j  §  /•  —  0 


UJ'P  — L O  o  —  0 


(30) 

(31) 

(32) 

(33) 


These  control  derivatives  will  ultimately  be  used  to  design  the  control  allocation  portion  of 
a  6  DOF  control  law  for  the  flapping  wing  MAY. 


VII.  Definite  Integrals  for  Computation  of  Time  Averaged 

Aerodynamic  Forces  and  Moments 

In  order  to  calculate  the  cycle-averaged  forces  and  moments,  it  will  be  necessary  to 
evaluate  numerous  integrals.  Many  have  no  indefinite  integral  solutions.  For  example,  many 
of  the  integrands  will  be  of  the  form  cos(costot)  or  sin  (cos  cot) .  These  terms  arise  because  of 
the  coordinate  transformations  of  forces  from  local  frames  to  the  body  frame  and  because 
of  the  fact  that  a  cosine  function  was  chosen  to  drive  the  position  of  the  wing.  Fortunately, 
definite  integrals  involving  such  functions  exist  over  the  intervals  of  interest  for  the  present 
problem  and  can  be  derived  from  results  presented  in  Gradshteyn  &  Ryzhik.6  The  solution 
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of  many  of  these  definite  integrals  involve  a  Bessel  function  of  the  first  kind,  as  well 

as  Struve  functions.  The  integrals  involving  Struve  functions  evaluate  to  zero  over  the  limits 
of  integration  that  arise  in  the  formulation  of  this  particular  problem.  For  convenience,  the 
solutions  to  the  definite  integrals  that  are  required  to  compute  the  cycle-averaged  forces  and 
moments  are  provided  below. 

h  =  I  sin2  [(a;  -  5)t]  dt  =  —  (34) 

Jo  2(cu  o) 


1-2  =  sin2  [(a;  —  8)t]  sin  (cos  [(a;  —  8)t])  dt  =  0 


h=  sin2  [(a;  —  5)t]  cos  (cos  [(a;  —  <5)f])  dt  = 


{u-S) 


h=  sin2  [(a;  +  a)t  +  £]  dt  = 


7 r 


2(u  +  a) 


J.5  =  /  sin2  [(a;  +  a)t  +  sin  (cos  [(a;  +  cr)t  +  ^])  dt  =  0 


Ie  =  /  sin2  [(a;  +  a)t  +  £]  cos  (cos  [(a;  +  cr)t  +  ^])  dt  — 


7rJi(l) 

u  +  a 


(35) 

(36) 

(37) 

(38) 

(39) 


VIII.  Cycle  Averaged  Forces 

In  this  section,  the  cycle-averaged  aerodynamic  forces  are  evaluated.  These  forces  are 
integrated  with  respect  to  time  as  each  wing  operates  over  an  upstroke  and  downstroke. 
Later,  the  forces  are  used  in  the  development  of  a  control  law  by  evaluating  the  change  in 
cycle-averaged  forces  with  respect  to  a  change  in  the  control  input  variables. 


A.  X  Force 


Substituting  the  expression  for  the  instantaneous  x-body  force  from  Table  3  into  Equation  29 
produces 

Urw 


vb  _ 
— 


2vr 


xJRW~°RW 


JRWU 


(■ t)dt  + 


JRWD 


(t)dt 


u’RW~SRW 


(40) 


Approved  for  public  release;  distribution  unlimited. 
17  of  25 


Substituting  Equations  6,  26,  and  27  into  Equation  40  yields 


X 


■B  kLOJRW 


RW 


2vr 

2? r 


ujRW~6RW 


°RW 


(urw  ~  $rw)2  sin2[(u;RW  —  $Rw)t]dt+ 
'o 

(ujrw  +  crw )2  sin2[(cdfw  +  aRw)t  +  I 


°RW  —°RW 


(41) 


Noting  that  the  time  varying  functions  under  the  integral  signs  are  of  the  form  R  and  I4  as 
given  by  Equations  34  and  37,  respectively,  Equation  41  can  be  written  as 


-^B  _  kjJ^jRW 

Xrw~^T 


(^rw  —  5rw)~Il  +  iyjRw  +  aRw)2h 


(42) 


or  simply 


X RW  —  L  m  (2 ujrw  —  &rw  +  ctrw) 


Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 

—B  k^LW  , 


X  lw  — 


-(2ulW  —  $LW  +  CTlw) 


(43) 


(44) 


Note  that  both  XRW  and  X LW  are  positive  quantities.  At  a  hover  condition,  where  the 
x-body  axis  is  normal  to  the  surface  of  the  earth,  the  forces  produced  by  both  wings  act 
to  counter  the  vehicle  weight.  It  is  also  worthwhile  to  note  that  because  of  the  geometric 
constraints  imposed  by  the  passive  wing  rotation  joints  shown  in  Figure  2,  the  instantaneous 
x-body  force  is  always  positive  or  zero  which  directs  the  force  along  the  positive  longitudinal 
axis  of  the  vehicle. 


B.  Y  Force 


In  order  to  compute  the  cycle-averaged  y-body  force  for  each  wing,  substitute  the  expression 
for  the  instantaneous  y-body  aerodynamic  force  from  Table  3  into  Equation  29  to  obtain 


—b  _  uRW 
Y  rw  ~ 


2vr 


VRW~°RW 


— DRWu(t )  sin [4>nw{t)]dt  + 


DRWD(t)  sin [<f>Rw(t)]dt 


°RW  —°RW 


5) 
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Substituting  Equations  6,  21,  22,  26,  and  27,  into  Equation  45  produces 


{urw  —  $rw)2  sin2[(u;RW  —  $Rw)t]  sin  (cos[(curw  —  $Rw)t])  dt+ 


(vrw  +  orw )2  sin2  1(lorW  +  <jRw)t  +  Crw]  sin  (cos[(u;Rvv  +  &Rw)t  +  Crw})  dt 

RW 

(46) 

Noting  that  the  time  varying  functions  under  the  integral  signs  are  of  the  form  /2  and  /5  as 
given  by  Equations  35  and  38  respectively,  Equation  46  becomes 

Y  rw  =  0  (47) 

Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 

YBlw  =  0  (48) 

Physically,  this  is  a  result  of  the  fact  that  the  component  of  drag  pointing  in  the  y-body 
direction  points  the  same  amount  of  time  in  the  positive  direction  and  at  the  same  magnitude 
as  it  points  in  the  negative  y-body  direction  over  the  course  of  a  wing  beat  cycle.  Note  that  if 
a  non-zero  cycle-averaged  force  in  the  y  direction  were  required,  the  frequency  of  the  driving 
oscillator  would  have  to  be  varied  over  finer  intervals  within  each  wing  beat  cycle.  In  other 
words,  the  wing  beat  frequency  would  have  to  be  changed  more  than  once  per  cycle. 


Y 


B 

RW 


knujRw 
2t r 


C.  Z  Force 


The  calculation  of  the  cycle-averaged  force  in  the  z-body  direction  follows  a  similar  procedure. 
Substituting  the  expression  for  the  instantaneous  z-body  aerodynamic  force  from  Table  3  into 
Equation  29  we  obtain: 


—B  _  UJrW 

z rw  — 


2tt 


wRW^5RW 


DRWuif)  cos [4>Rw(t)]dt  + 


■ DRWD{t )  cos [(f>Rw(t)]dt 


WRW~SRW 


9) 
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Substituting  Equations  6,  21,  22,  26,  and  27  into  Equation  49  yields 


'  RW 


ko^RW 
27 r 

2t r 


dRW~°RW 


VRW 


dRW  —°RW 


{urw  ~  Srw)2  sin2[(u;Rw  —  $Rw)t\  cos  (cos[(curvv  —  <5 Rw)t ])  dt— 

'  0 

(wrvk  +  aRw)2  sin2 [(ogRve  +  ^ Rw)t  +  £rw]  cos  (cos[(curvv  +  crRW)t  +  Crw})  dt 

(50) 


Noting  that  the  time- varying  functions  under  the  integral  signs  are  of  the  form  J3  and  J6  as 
given  by  Equations  36  and  39  respectively,  Equation  50  can  be  written  as 


7B 
^ RW 


kpUjRw 

2n 


(urw  —  5rw)2Iz  ~  (urw  +  &rw)2I& 


(51) 


or  simply 

(52) 

Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 

zlw  =  =^m^L{SLW  +  aLW)  (53) 

The  result  for  the  z-body  axis  is  important  because  it  means  that  nonzero  cycle-averaged 
forces  can  indeed  be  generated  and  used  to  induce  fore  and  aft  linear  accelerations.  It  can 
also  be  used  to  generate  rolling  moments. 


IX.  Cycle  Averaged  Moments 

In  this  section,  expressions  for  the  cycle-averaged  aerodynamic  moments  are  derived.  The 
cycle-averaged  moments  are  later  used  in  the  development  of  a  control  law  that  requires  an 
evaluation  of  the  change  in  cycle-averaged  moments  with  respect  to  a  change  in  the  control 
input  variables. 


A.  Rolling  Moment 

Substituting  the  expression  for  the  instantaneous  x-body  moment  from  Equations  15  and  16 
into  Equation  29  gives 


-B 


Mxrw  ~ 


Urw 

2tt 


wRW~SRW 


M- 


xRWU 


(t)dt  + 


M. 


xRWD 


(t)dt 


WRW~SRW 


(54) 
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Substituting  Equations  6,  21,  22,  26,  and  27  into  Equation  54  yields: 


wB  -  URW 

1V1XRW  ~ 


2vr 


w 


kDy%,P{uRW  ~  SrwY 


VRW~°RW 


k-D^^RW  ~  Srw)2 
ko^ZftiuRw  —  Srw 


VRW  —°RW 


°RW  —°RW 


sin2[(otmy  —  dRw)t]dt  +  •  •  • 

Jo 

cos  {cos  [(c Jrw  ~  <5rvu)£]}  sin2[(uiRW  —  SRw)t]dt  +  . . . 
sin  {cos  [( urw  —  SRw)t}}  sin2[(u;,Rn/  —  8Rw)t]dt  —  . . . 


^d^rw  +  ORwfy™/  /  sin2  [(ujrw  +  crRw)t  +  Crw]  dt  —  . . . 


UJRW~SRW 
2-7 r 


koi^Rw  +  (JrwY 


sin2  [{urw  +  o' R\y)t  +  Crw] 


°RW  —°RW 


cos  {cos  {(u)rw  +  orw)^  +  £rvk]}  dt  —  . . . 

2-7 r 

ko(^RW  +  <Jrw)2Azr  f  sin2  [(urw  +  &Rw)t  +  £rw] 


VRW  —°RW 


sin  {cos  [( urw  +  <JRw)t  +  £rw}}  dt] 


(55) 


Note  that  the  integrals  are  of  the  form  given  in  Equations  34  -  39;  thus,  Equation  55  can  be 
written  as 


-jj-b  _  ^Rwk'D 
1V1xRW  ~ 


2vr 


- 5 


RW 


vZPh  +  A  4/2  +  ^I: 


IU 


(urw  +  ctrw )2  ylPh  +  A^/,5  +  — /e  | 

Substituting  the  results  for  the  definite  integrals,  I\  -  /6,  and  simplifying  yields 

^D^RW  (8r\V  +  &Rw) 


TTb 
MXRW  = 


[yZP  +  wM  1)] 


Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 

-ir-j-B  kR)UJLw  (8lw  +  &lw) 


Mxlw  ~ 


[yZP  +  wJ i(!)] 


(56) 


(57) 


(58) 


Note  that  without  split-cycle  frequency  modulation,  i.e.,  5rw  =  ctrw  =  8lw  =  &lw  =  0,  it 
would  not  be  possible  to  generate  non- zero  cycle-averaged  rolling  moments  on  this  aircraft. 
Thus,  because  of  the  constraints  resulting  from  the  minimally  actuated  vehicle  design,  split- 
cycle  modulation  is  a  sufficient  condition  for  generating  rolling  control  moments  for  the 
proposed  aircraft. 
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B.  Pitching  Moment 

The  instantaneous  y-body  moment  is  determined  by  substituting  Equations  15  and  16  into 
Equation  29 


WB  =^RW 
IWyRW  2n 


dRW~°RW 


MyRwu(t)dt  + 


MyRwoifydt 


'0 


dRW  —°RW 


(59) 


Substituting  Equations  6,  21,  22,  26,  and  27  into  Equation  59  yields: 

_ 7 T _ 

RW~SRW  .  2 


WB  =^RW 
1  yRW  2tt 


-kD  AxR 


‘1RW~°RW 


-kj/x^p  cos  a 


r  WP  ■ 
—koxr.n  sm  a 


{urw  —  $rw)2  kRAzR  I  siii2 [(u> rw  —  SRw)t]dt  +  ... 

Jo 

sin2  [(a jrw  —  SRw)t]  cos  {cos(curw  —  $Rw)t}  dt  +  . . . 

sin2  [(ujrw  ~  $Rw)t]  cos  {cos  [(a jRw  ~  d Rw)t ]}  dt  +  . . . 


°RW~°RW 


UJRW~SRW 


sin2[(agRw  —  SRw)t]  cos  {cos  [{ujrw  —  SRw)t]}  dt  +  . . . 


kLy, 


WP 

cp 


dRW  —°RW 


sin2[(a;.RVF  —  $Rw)t]  sin  {cos  {(ivrw  —  dRw)t]}  d^ 


+ 


(ijjRW  +  CrwY 


kLAzp 


sin"  [(ujrw  +  &Rw)t  +  £rw]  dt  +  . . . 


" RW~SRW 


kpAx 


sin2  [(a )rw  +  crRw)t  +  £rw] 


u1RW-5RW 

cos  {cos  [(uRW  +  (jRw)t  +  £,Rw]}dt  +  ... 


kLxZP  cos  a 


sin2  \(uRW  +  crRw)t  +  £rvc] 


°RW~°RW 


COS  {cos  [(a Jrw  +  o Rw)t  +  £rw}}  dt  +  . . . 

2?r 

sin"  \{uRW  +  crRw)t  +  £rw] 


kpxYp  sin  a 


dRW~°RW 


cos  {cos  [(uRW  +  aRW)t  +  £Rw\}dt  +  ... 

2ir 

sin2  [{ujrw  +  &Rw)t  +  £rw] 


kiycv 


WP 


°RW  —°RW 


sin  {cos  [U>;ic  +  aRW)t  +  £rw]}  dt]} 


(60) 
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Note  that  the  integrals  are  of  the  form  given  in  Equations  34  -  39;  thus,  Equation  60  can  be 
written  as 


Mvrw  -  <W)2  [kLAz%Ii  -  koAx^ls  -  kLx^pp  cos  al3— 

kDx™p  sin  al3  +  kLy^pI2]  +  (uRW  +  aRW )2  [kLAz%I4+  (61) 

kD AxRI6  +  kLx^p  cos  a/6  +  kDx™p  sin  al6  +  kLy^pI5\  } 


Substituting  the  results  for  the  definite  integrals,  I\  -  h,  and  simplifying  produces 


-jj-B  _  kLX>RW 

mvrw  ~  2 


Azp 


xcP  ^i(l)  coso;(5rw  +  &Rw)  4 - t^-(2urw  +  <xrw  ~  Srw) 


koJiiXt^Rw 


(62) 


{  (Srw  +  cfrw)  [x™P  sin  a  +  AxR]  } 


Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 


-j-j-B  _  kRUJLW 

My  lw  ~  2 


z  AzB 

x^P  J i(l)  cos a(dLW  +  <jlw)  4 - y^(2 ulw  +  &LW  ~  $lw) 


+ 


koJliXfX’LW 


(63) 


{ (hw  +  ^ir)  [x™pP  sin  a  +  Axf  ]  } 


Note  that  without  split-cycle  frequency  modulation,  i.e.,  5rw  —  &rw  —  &lw  —  &lw  —  0, 
there  still  exists  a  non-zero  cycle-averaged  pitching  moment  if  A zB  ^  0  and  A zB  ^  0. 
This  result  yields  further  insight  into  vehicle  design  for  controllability.  It  suggests  that  the 
wing  root  hinge  point  should  be  placed  such  that  its  z-body  location  is  coincident  with  the 
nominal  vehicle  center-of-gravity,  i.e.,  A zR  =  0,  Azf  =  0.  Such  an  arrangement  will  yield 
a  zero  cycle-averaged  pitching  moment  when  the  split-cycle  parameters  are  zero,  which  is 
a  desirable  feature  for  maintaining  hover.  Note  that  the  current  Harvard  RoboFly  aircraft 
does  not  incorporate  this  feature  and  would  immediately  pitch  forward  if  not  constrained 
by  the  guide  wires  used  in  their  experiment.  Later,  the  effect  of  a  bob-weight  actuator,  that 
can  be  used  to  manipulate  pitching  moments,  will  be  examined. 


C.  Yawing  Moment 

The  instantaneous  z-body  moment  is  determined  by  substituting  Equations  15  and  16  into 
Equation  29  to  obtain 


Mzrw  ~ 


Urw 

2tt 


MzRwu(t)dt  + 


Mz  RWD(t)dt 


(64) 
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Substituting  Equations  5,  21,  22,  26,  and  27,  into  Equation  64  yields: 


Mzrw  =  \  -  $rw)2  ~kL^  /  RW  RW  sin2 [(coRW  ~  SRW)t\dt 


kD 


sin2[(cuRvv  —  SRw)t]  sin  {cos(curvv  —  SrwY}  dt  —  .  . . 


WP  [  ^RW~5RW  .  2  r/  r  \  t  r  r/  c-  \n7 

kLycp  /  snr[(ww  -  0RW)t\  cos  {cos  [(cuw  -  dRW)t\ }  dt  -  ... 

Jo 


kLx'jp  cos  a 


sin2[(<^Riy  —  dRw)t\  sin  {cos  [(ujrw  ~  SrwY]}  dt  —  .. . 


kDx^p  sin  a  /  ™  sin2[(o>w  -  dRW)t]  sin  {cos  [(c vRW  -  5RW)t]}dt  +  ... 


(ujRW  +  (TrwY  ~kL~2 


2  ,  W  I  “RW  ■  2 


sin  [{ojrw  +  <JRw)t  +  t,Rw\  dt  +  . . . 


kDAxft 


UJRW~°RW 

2? r 

WP  [“RW 

cp  I 

w RW~SRW 


sin2  [{wrw  +  &Rw)t  +  £rw]  sin  {cos  [(ujrw  +  ctrwY  +  £rw]}  dt  —  . 


sin2  [( Urw  +  crRW)t  +  £rw]  cos  {cos  [{(jJrw  +  <JRw)t  +  £rw}}  dt  + 


cos  a  /  sin"  [{urw  +  aRw)t  +  £mv] 


sin  {cos  [(c urw  +  <JRW)t  +  Crw}}  dt  + 


koxZ/  sin  a 


sin2  [(a )rw  +  &Rw)t  +  £rw] 


WRW~SRW 


sin  {cos  [{uRW  +  aRW)t  +  Crw}}  dt}} 


Note  that  the  integrals  are  of  the  form  given  in  Equations  34-39;  thus,  Equation  65  can  be 


written  as 


B  OJrw 


MZ  DTI/  - 


-kL—h  -  kDAxpI2  -  kLy™pI3- 


r  iv 

kLxYpP  cos  al2  -  kDxYpP  sin  al2]  +  (uRW  +  crRW)2  ~kL—I4+ 
kDAxpI5  -  kLyYp  Ph  +  kLx™pP  cos  a/5  +  kDxYpP  sin  al5]  } 


Substituting  the  results  for  the  definite  integrals  Ii-Iq,  and  simplifying  yields 


_—kLOJRW 
MzRW  —  /> 


](yZPM  i) 


+  (2 ujrw  —  5rw  +  &rw) 
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Following  a  similar  procedure  for  the  left  wing,  it  can  be  shown  that 


■j~rB  _  ^'LULW 

1V1zLW  ~  o 


(yWpPJ i(l)  +  —  j  (2a >RW  —  5rw  +  ctrw) 


(68) 


Without  split-cycle  frequency  modulation,  there  exists  a  non- zero  cycle-averaged  yawing 
moment  on  each  wing  since  ylpp  ^  0  and  w  ^  0.  When  the  split-cycle  parameters  for  each 
wing  are  zero,  the  cycle-averaged  moments  are  opposing  and  balance  one  another.  Since 
the  fundamental  wing  beat  frequency,  iv,  can  be  independently  varied  for  each  wing,  yawing 
moments  can  be  generated  without  varying  the  split-cycle  parameters. 


X.  Summary 

The  analysis  presented  shows  that  the  use  of  split-cycle  constant-period  frequency  modu¬ 
lation,  with  independently  actuated  wings,  allows  manipulation  the  x-body  and  z-body  axis 
forces.  The  y-body  force  is  not  directly  controllable  using  the  split-cycle  approach  as  pre¬ 
sented.  Direct  manipulation  of  the  y-body  force,  given  the  physical  design  constraints  of  the 
proposed  vehicle,  would  require  that  the  frequency  of  the  oscillators  driving  the  wings  be  var¬ 
ied  more  than  once  per  cycle.  Rolling  and  yawing  moments  can  also  be  independently  varied 
using  the  split-cycle  technique.  A  pitching  moment  can  be  generated;  however,  it  cannot  be 
generated  independently  of  the  z-body  force  or  the  rolling  moment.  In  Part  II,  a  bob-weight 
is  introduced  to  cancel  these  undesirable  changes  in  the  pitching  moment  when  manipulating 
the  split-cycle  parameters  to  generate  z-body  forces  and  rolling  moments.  Given  the  ability 
to  independently  manipulate  5  out  of  the  6  cycle-averaged  body-axis  forces  and  moments, 
untethered  controlled  flight  with  insect-like  maneuverability  appears  to  be  feasible. 
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